/********************************************************************************
do-file to replicate tables in:
	Chacon, Mario and Jeff Jensen. 2019. "Democratization, De Facto Power, and Taxation: 
	Evidence from Military Occupation during Reconstruction" World Politics 

data file needed: Dataset_1 (Main Panel) 
*/

clear all 
set more off

**Set Working Directory 
*cd "C:\****\ Replication Files (Chacon & Jensen, WP 2019)\Datasets\"

/*
***	Plotted estimates coming from the following regressions: 
use Dataset_1, clear 

   foreach num of numlist 1870 1880 1890 {
    reg ln_state_tax_real_pc black_sh_1860 lnTotPop1860 lnFarValRPc1860 if sample==1 & year==`num' & d_occup_recon==0, beta
	reg ln_state_tax_real_pc black_sh_1860 lnTotPop1860 lnFarValRPc1860 if sample==1 & year==`num' & d_occup_recon==1, beta 
	    }

foreach num of numlist 1870 1880 1890 {
    reg ln_county_tax_real_pc black_sh_1860 lnTotPop1860 lnFarValRPc1860 if sample==1 & year==`num' & d_occup_recon==0, beta
	reg ln_county_tax_real_pc black_sh_1860 lnTotPop1860 lnFarValRPc1860 if sample==1 & year==`num' & d_occup_recon==1, beta 
	    }		
*/

use Dataset_Fig_3, clear

twoway (line state_cont year, lcolor(black) lpattern(dash) subtitle("ln(Per capita State Tax Revenue)") ytitle("")) ///
		(scatter state_cont year, msymbol(oh) msize(vlarge) mcolor(black)) ///
       (line state_treat year, lcolor(black)) ///
	   (scatter state_treat year, msymbol(o) msize(vlarge) mcolor(black)), ///
	   name(A, replace) scheme(s1mono) legend(order(2 4) region(fcolor(white) lcolor(none))) ///
	   xline(1877, lpattern(tight_dot)) 

twoway (line county_cont year, lcolor(black) lpattern(dash) subtitle("ln(Per capita County Tax Revenue)") ytitle("")) ///
		(scatter county_cont year, msymbol(oh) msize(vlarge) mcolor(black)) ///
       (line county_treat year, lcolor(black)) ///
	   (scatter county_treat year, msymbol(o) msize(vlarge) mcolor(black)), /// 
	   name(B, replace) scheme(s1mono) legend(order(2 4) region(fcolor(white) lcolor(none))) ///
	   xline(1877, lpattern(tight_dot)) 
		
graph combine A B, scheme(s1mono) 

